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The latest determination of the extragalactic gamma-ray background (EGRB) radiation by Fermi 
is compared with the theoretical prediction of the blazar component by Inoue & Totani (2009; here- 
after IT09). The Fermi EGRB spectrum is in excellent agreement with IT09, indicating that blazars 
are the dominant component of the EGRB, and contributions from any other sources (e.g., dark 
matter annihilations) are minor. It also indicates that the blazar SED (spectral energy distribu- 
tion) sequence taken into account in IT09 is a valid description of mean blazar SEDs. The possible 
contribution of MeV blazars to the EGRB in the MeV band is also discussed. In five total years 
of observations, we predict that Fermi will detect ~1200 blazars all sky down to the correspond- 
ing sensitivity limit. We also address the detectability of the highest-redshift blazars. Updating 
our model with regard to high-redshift evolution based on SDSS quasar data, we show that Fermi 
may find some blazars up to z ~ 6 during the five-year survey. Such blazars could provide a new 
probe of early star and galaxy formation through GeV spectral attenuation signatures induced by 
high-redshift UV background radiation. 



INTRODUCTION 



The origin of the extragalactic diffuse gamma-ray 
background (EGRB) is a long-standing puzzle. First 
discovered by the SAS 2 satellite [l|, y], its exis- 
tence was subsequently confirmed up to ~ 50 GeV 
by EGRET (Energetic Gamma-Ray Experiment Tele- 
scope) on board the Compton Gamma Ray Obser- 
vatory. EGRET measured a flux of about lx 10~ 5 
photons cm~ 2 s _1 sr" 1 above 100 MeV and an ap- 
proximately power-law spectrum with photon index of 

2 over the wide range of 30 MeV - 50 GeV @, i]. 

Very recently, Ref. [f| reported Fermi observations of 
the EGRB spectrum up to 100 GeV, which connects 
smoothly to the EGRET results below - 200 MeV 
0| . However, above this energy, the spectrum is dis- 
crepant from EGRET, being a single power-law with 
index ~ —2.4, and showing no signs of a GeV excess. 

Although different types of gamma-ray sources 
(e.g., clusters of galaxies or dark matter annihila- 
tion) have been proposed to be significant contrib- 
utors to the EGRB [see, e.g., Ref. @] and refer- 
ences therein], active galactic nuclei (AGNs) of the 
blazar class are considered the primary candidates, 



since almost all of the extragalactic sources detected 
by EGRET were blazars. The blazar contribution to 
the EGRB has been estimated by numerous authors, 
who have reached different conclusions using different 
approaches, ranging from 20% to 100 % of the ob- 
served flux [6h18|. 

In most past studies, the spectral energy distri- 
butions (SEDs) of the blazars were assumed to be 
power laws for all objects. In such models, it is ob- 
vious that the predicted EGRB spectrum is mainly 
determined by the assumed power-law indices. How- 
ever, multi-wavelength observational studies indicate 
the existence of a non-trivial trend among blazar 
SEDs: the energies of the two characteristic spectral 
peaks in blazar SEDs (each likely reflecting the syn- 
chrotron and inverse Compton emission) systemati- 
cally decrease as the bolometric luminosity increases 
|19[-[25|. This is often referred to as the blazar SED 
sequence. Although its validityis currently still a mat- 
ter of debate (e.g., Ref. [2~4l. l25l. l27l|). one can make 
a non-trivial prediction of the EGRB spectrum if this 
blazar sequence is assumed. 

Here we compare the latest Fermi EGRB results 
with our model predictions, previously published in 
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Ref. [26| (IT09) before the former was announced. 
IT09 calculated the EGRB from blazars by construct- 
ing a blazar gamma-ray luminosity function (GLF) 
model that is consistent with the flux and redshift dis- 
tributions of the EGRET blazars, accounting for the 
blazar sequence as well as the luminosity-dependent 
density evolution (LDDE) scheme that describes well 
the evolution of the X-ray luminosity function (XLF) 
of AGNs. By introducing the blazar SED sequence, 
we were able to make a reasonable and non-trivial pre- 
diction of the EGRB spectrum that can be compared 
with observations including the new Fermi data. 

Recently, Ref. [H has shown that the EGRB in 
the MeV band can be naturally explained by normal 
(i.e., non-blazar) AGNs that compose the cosmic X- 
ray background. If so, they should also contribute to 
the EGRB at < 1 GeV. We will investigate what frac- 
tion of the observed EGRB can be explained by the 
sum of the blazar and non-blazar AGN components. 
The contribution of blazars to the EGRB in the MeV 
band is also discussed and compared with a recent 
study on the possible relevance of the so-called "MeV 
blazars" [H. 

We will then make quantitative predictions for 
the Fermi Gamma-ray Space Telescope [30l | for its 
five-year survey. The first Fermi catalog for bright 
gamma-ray sources including AGNs has recently been 
released [3l], [H[ . The number of currently detected 
blazars is about 100, larger than that of EGRET by 
a modest factor of about 2. We limit the discussion 
of the blazar data to only the EGRET sample in this 
work, to be considered as a theoretical prediction from 
the pre-Fermi era. In the future, a much larger num- 
ber of blazars should be detected with Fermi, which 
can be compared further with our predictions. 

We also discuss the highest-redshift blazars de- 
tectable by Fermi, by updating the IT09 GLF model 
with respect to high-redshift evolution. We briefly 
mention the possibility of utilizing such blazars to 
obtain information on the high-redshift extragalactic 
background light and early galaxy formation through 
spectral absorption features. 

Throughout this paper, we adopt the standard cos- 
mological parameters of (h, £Im, fiA)=(0.7,0.3,0.7). 

II. THE BLAZAR SED SEQUENCE 

Ref. 0, 0, [23[ constructed an empirical blazar 
SED model to describe the SED sequence, based on 
fittings to observed SEDs from radio to 7-ray bands. 
These models are comprised of the two components 
(synchrotron and IC), and each of the two is de- 
scribed by a linear curve at low photon energies and 
a parabolic curve at high energies. 

We construct our own SED sequence model mainly 
based on the SED model [H[ , because there is a math- 
ematical discontinuity in the original model of Ref. 
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FIG. 1: The blazar SED sequence. The data points 
are the average SED of the blazars studied by Ref. 
[20L [2^ |. The solid curves are the empirical SED se- 
quence models constructed and used in this paper. The 
model curves corresponds to the bolometric luminosities 
of log 10 (P/erg s" 1 ) = 49.50, 48.64, 47.67, 46.37, and 45.99 
(from top to bottom). 

[25 | . Our own SED sequence formula is described in 
the Appendix of Ref. [26( in detail. In FigdJ we show 
this empirical blazar SED sequence model in compar- 
ison with the observed SED data [10,111]. 



III. THE MODEL OF GAMMA-RAY 
LUMINOSITY FUNCTION OF BLAZARS 

The cosmological evolution of AGN XLF has been 
investigated intensively [33j-|36j. These studies re- 
vealed that AGN XLF is well described by the LDDE 
model, in which peak redshift of density evolution 
increases with AGN luminosity. Here we construct 
blazar GLF models based on the two XLFs derived 
by Ref. [33[ (hereafter U03; in hard X-ray band) and 
Ref. [34| (hereafter H05; in soft X-ray band). The use 
of LDDE in blazar GLF has been supported from the 
EGRET blazar data [|. 

We simply assume that the bolometric luminosity 
of radiation from jet, P, is proportional to disk X-ray 
luminosity, Lx ■ We relate these two by the parameter 
q, as P = 10 q Lx- Here, we define the disk luminosity 
Lx to be that in the rest-frame 2-10 and 0.5-2 keV 
bands for the U03 XLF and the H05 XLF , respec- 
tively. Thus, the luminosity at rest 100 MeV, L 7 , and 
Lx have been related through P. 

The blazar GLF p 1 is then obtained from the AGN 
XLF, p x , as p 7 (Lj,z) = K^-p x (L x ,z), where p 1 
and px are the comoving number densities per unit 
gamma-ray and X-ray luminosity, respectively. The 
parameter k is a normalization factor, representing 
the fraction of AGNs observed as blazars. See the §3 
in Ref. [26j for a detail. 

We use the maximum likelihood method to search 
for the best-fit model parameters of the blazar GLF 
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to the distributions of the observed quantities of the 
EGRET blazars (gamma-ray flux and redshift). See 
Ref. 26] for a detail. We take q and the faint-end 
slope index of XLF 71 in addition to q as free param- 
eters . These are hereafter called as U03(g), H05(q), 
U03(<7, 71) and H05(g, 71) fits, respectively. Figure 
[5] shows the distributions of redshift and gamma-ray 
luminosity predicted by the best-fit models, in com- 
parison with the EGRET data. 

We performed the Kolmogorov-Smirnov (KS) test 
to see the goodness of fits for the best-fit results of 
each of the four fits, and the chance probabilities of 
getting the observed KS deviation. Since the best KS 
test value is obtained for the U03(g,7i) GLF model, 
we use this GLF model as the baseline below. 



IV. THE EGRB SPECTRUM 

We calculate the EGRB spectrum by integrating 
our blazar SED sequence model using the blazar GLF 
derived in mill Since EGRET has already resolved 
bright gamma-ray sources, we should not include 
those sources in the EGRB calculation. Therefore, 
the maximum gamma-ray luminosity in the integra- 
tion should be i 7 (i 7 EGRET, z), which is the luminosity 
of the blazar having a flux Ky = -Fegret at a given 
z, where F E g RET = 7x 1CP 8 photons cm~ 2 s _1 is the 
EGRET sensitivity above 100 MeV. 

High energy photons (> 20 GeV) from high redshift 
are absorbed by the interaction with the cosmic in- 
frared background (GIB) radiation (37l - l40| . We adopt 
the model of Ref. [38] for the CIB optical depth. We 
also take into account the cascading emission to calcu- 
late EGRB spectrum, considering only the first gener- 
ation of created pairs [13] • We found that the amount 
of energy flux absorbed and reprocessed in intergalac- 
tic medium is only a small fraction of the total EGRB 
energy flux , and hence the model dependence of CIB 
or the treatment of the cascading component does not 
have serious effects on our conclusions in this work 

In addition to blazars, we take into account non- 
blazar AGNs as the source of EGRB. Ref. [II] (here- 
after ITU08) has shown that non-blazar AGNs are 
a promising source of EGRB at 10 MeV. In this 
scenario, a nonthermal power-law component extends 
from hard X-ray to —10 MeV band in AGN spec- 
tra, because of nonthermal electrons that is assumed 
to exist in hot coronae around AGN accretion disks. 
The existence of such nonthermal electrons is theo- 
retically reasonable, because magnetic reconnection is 
the promising candidate for the heating source of the 
hot coronae [461 ] where hard X-ray photons are emit- 
ted. Magnetic reconnections should produce nonther- 
mal particles, as is well known in e.g., solar flares or 
the earth magnetosphere. Although several sources 
have been proposed as the origin of the MeV back- 



ground, this model gives the most natural explanation 
for the observed MeV background spectrum that is a 
simple power- law smoothly connected to CXB. 

Theoretically, there is no particular reason to ex- 
pect a cut-off of the nonthermal emission around 10 
MeV, and it is well possible that this emission ex- 
tends beyond 10 MeV with the same power index. 
Then, this component could make some contribution 
to GeV EGRB. The EGRB spectra from the non- 
blazar AGNs are calculated based on the model of 
ITU08. The model parameters of ITU08 have been 
determined to explain the EGRB in the MeV band. 
There is some discrepancy between the MeV EGRB 
data of SMM and COMPTEL [44], and the rea- 
son for this is not clear. Here we use two model pa- 
rameter sets of (r,7 tr ) = (3.5,4.4) and (3.8,4.4) in 
ITU08, where T is the power-law index of nonthermal 
electron energy distribution and 7 tr is the transition 
electron Lorentz factor above which the nonthermal 
component becomes dominant. These two parameter 
sets are chosen so that they fit to the COMPTEL and 
SMM data, respectively. 



A. On the Origin of the EGRB: Comparison 
with the Fermi EGRB Spectrum 

Fig. [3] compares our model predictions with the ob- 
served EGRB data, including the very recent Fermi 
results |[ [Hi]. Since the flux thresholds of EGRET 
and Fermi are different, we must be careful when 
directly comparing with the Fermi data. However, 
when the sensitivity threshold of Fermi is set to be 
the same as that of EGRET, the EGRB flux changes 
by only -10 % [13]. Therefore, our predicted EGRB 
spectrum can be considered to be in excellent agree- 
ment with the Fermi data. This agreement strongly 
suggests that blazars are the principal sources of the 
EGRB above 1 GeV, and the contribution from other 
sources like dark matter annihilation is minor. This 
also implies that use of the blazar sequence is a valid 
assumption. 

The EGRB prediction using the ITU08 non-blazar 
background model with T = 3.5 also agrees nicely 
with the observed data from X-rays to 100 GeV. In 
this case, the predicted EGRB flux above 100 MeV 
can account for 80% of the observed flux, consider- 
ably higher than in previous studies [1, H, fl3j ] . It 
should be noted that the contribution to the EGRB 
flux above 100 MeV from blazars alone is ~45 %, while 
the remaining ~35 % comes from non-blazar AGNs. 

Based on these results and considerations, we con- 
clude that most, and probably all, of the EGRB flux 
can be explained by blazars together with non-blazar 
AGNs, whose luminosity functions are consistent with 
the EGRET blazar data and X-ray AGN surveys. 
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FIG. 2: Redshift and gamma-ray luminosity (yL v at 100 MeV) distributions of EGRET blazars. The histogram is the 
binned EGRET data, with one sigma Poisson errors. The four model curves are the best-fits for the GLF models of 
U03(g), H05(g), U(g,7i), and H05(g, 71), as indicated in the figure. Bottom panels: the same as the top panels, but for 
cumulative distributions. 
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FIG. 3: The EGRB spectrum from non-blazar AGNs and 
blazars for the GLF model of U03(q, 71). The model curves 
for the blazar component (absorbed+cascade), non-blazar 
AGN component, and the sum of the two are shown. Note 
that two models are plotted for the non-blazar component 
with different values of Y. The observed data of HEAO- 
1 E3, Swift/BAT gH, SMM @, COMPTEL [p and 
EGRET 0, [5[ are shown. We also plot an alternative ver- 
sion of the EGRET data denoted as "S98+S08" which is 
the original EGRET results modified with the correc- 
tion factors of Ref. (45|. The Fermi EGRB data is also 
shown [a|. 



B. Non-blazar AGNs and Blazars 

As can be seen in Fig. [31 the spectrum of the cos- 
mic MeV background connects smoothly with that of 
the cosmic X-ray background (CXB), while the spec- 
tral index of the EGRB changes abruptly around ~10 
MeV. This indicates that the sources of the EGRB 
below 10 MeV is the same population as that of the 
CXB (i.e. non-blazar AGNs), and the sources of the 
EGRB above 10 MeV is a different population (i.e. 
blazars). 



Recently, based on observations of Swift/BAT- 
selected blazars at 15-55 keV, Ref. [29[ suggested 
that MeV blazars contribute significantly to the MeV 
background. However, this conclusion is based on 
an extrapolation of blazar SEDs from 55 keV to the 
MeV band, and some fine tuning may be required in 
the MeV blazar SED to reproduce the EGRB spec- 
trum that connects smoothly with the CXB. More- 
over, the spectral break in the EGRB around 10 MeV 
implies that such MeV blazars are a distinct popula- 
tion from those contributing to the EGRB above 10 
MeV. Therefore, it seems more natural that the dom- 
inant contribution to the MeV background is coming 
from non-blazar AGNs as proposed by 28]. Further 
observational studies are desired to reveal the true ori- 
gin of the MeV background. 



V. PREDICTIONS FOR THE FERMI MISSION 



A. Expected Number of Blazars and Non-blazar 
AGNs 



The left panel of Figure 2] shows the cumulative dis- 
tribution of > 100 MeV photon flux of blazars, non- 
blazar AGNs, and the total of the two. The GLF 
predicts that about 750 and 1200 blazars should be 
detected by Fermi where we assumed the Fermi sensi- 
tivity to be Fii m = 3 x 10 -9 and 1 x 10 -9 photons cm~ 2 
s~ x above 100 MeV. These sensitivities correspond to 
~1 and ~5 years survey sensitivities. 

The model of Ref. [H|, and the PLE and LDDE 
models @] predicted ~10000, 5000, and 3000 blazars 
for the 1-year sensitivity limit, respectively. It is re- 
markable that the GLF models in this work predict 
significantly smaller numbers of blazars than most of 
previous studies. 

About 4-50 non-blazar AGNs would be detected 
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FIG. 4: Left: the cumulative flux distributions of blazars, non-blazar AGNs, and the total of the two. The detection 
limits of EGRET and Fermi are also shown. Right: the same as the left panel, but showing differential flux distribution 
multiplied by flux F ly to show the contribution to EGRB per logarithmic flux interval. The two models of non-blazar 
AGNs with different values of Y are shown. The total of blazar and non-blazar counts is also shown. 



by Fermi with 5-year survey, by their soft nonthermal 
emission from nonthermal electrons in coronae of ac- 
cretion disks, giving an interesting test for our model. 
We estimated the > 100 MeV gamma-ray flux from 
the observed hard X-ray flux of NGC 4151 (the bright- 
est Seyfert galaxy in all sky [35]) using the ITU08 
model [26(. Then, NGC 4151 is marginally detectable 
when r = 3.8, and easily detectable with T = 3.5 with 
one year survey sensitivity. 

The right panel of Figure H] shows the differential 
flux distribution of gamma-ray blazars and non-blazar 
AGNs multiplied by flux, showing the contribution to 
the EGRB per unit logarithmic flux interval. EGRB 
from blazars should practically be resolved into dis- 
crete sources. We find that the fraction of EGRB 
flux that should be resolved is 98% against the to- 
tal blazar EGRB flux in 5-year survey. On the other 
hand, EGRB from non-blazar AGNs is hardly resolved 
into discrete sources by Fermi. The predicted fraction 
of non-blazar EGRB flux resolved by Fermi is ^0.01 % 
with 5-year survey, respectively, against the total non- 
blazar EGRB flux. These results are in sharp contrast 
to those for blazars, although the contributions to the 
total EGRB by the two populations are comparable 
at around 100 MeV. 

The above results mean that a considerable part 
of EGRB at around 100 MeV will remain unresolved 
even with the Fermi sensitivity, because there is a 
considerable contribution from non-blazar AGNs to 
EGRB at ^100 MeV. However, the contribution from 
non-blazars should rapidly decrease with increasing 
photon energy and almost all of the total EGRB flux 
at > 1 GeV should be resolved into discrete blazars 
by Fermi, if there is no significant source contribut- 
ing to EGRB other than blazars. It should be noted 
that these predictions can be tested by Fermi rela- 
tively easily, without follow-up or cross check at other 
wavebands, once measurements of source counts and 
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FIG. 5: Expected cumulative redshift distribution of 
blazars detectable by Fermi above 100 MeV. Red and 
blue curves correspond to the flux sensitivity limits of i*Ii m 
(MOOMeV) = 3 x 10" 9 and 1 x 10" 9 in units of photons 
cm -2 s _1 . Solid and dashed curves correspond to whether 
or not SDSS quasar constraints are taken into account. 



EGRB flux have been done by the Fermi data. There- 
fore this gives a simple and clear test for our theoret- 
ical model |26[ . 



B. High Redshift Fermi Blazars 

Since Fermi is the most sensitive GeV gamma-ray 
observatory to date, it has the potential to discover 
many new high-redshift blazars. GeV photons from 
high-redshift sources are expected to suffer intergalac- 
tic absorption due to e± pair production interactions 
with the cosmic U V background radiation. The result- 
ing features in the spectra of such blazars could there- 
fore be a key probe of the poorly-understood UV back- 
ground radiation and provide valuable insight into the 
cosmic dark ages |5fj|-[5a|. 
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The high-redshift evolution of the GLF is uncer- 
tain, as there are no appropriate samples of gamma- 
ray blazars and X-ray AGNs above z — 3. Thus we 
constructed new GLFs by combining the U03 XLF 
with evolutionary constraints from the Sloan Digital 
Sky Survey (SDSS) quasar luminosity function data 
[531 ] . The SDSS quasar luminosity function is well 
constrained up to z ~ 5. 

Fig. [5] shows the expected cumulative redshift 
distribution of blazars that are detectable by Fermi 
above 100 MeV. Two sets of curves are plotted corre- 
sponding to the one- and five-year survey sensitivity. 
Solid and dashed curves refer to the cases of whether 
or not the SDSS quasar constraint is taken into ac- 
count. Note that intergalactic absorption effects are 
not included here, as they are not expected to be im- 
portant below 1 GeV 50^521] . The results imply that 
Fermi may detect some blazars up to z ~ 6 during the 
five-year survey. Measurements of the spectra of such 
blazars above 1 GeV through subsequent, deep obser- 
vations should provide interesting constraints on the 
high-z UV background. More details will be presented 
in a forthcoming publication. 



VI. CONCLUSIONS 

We have compared our prediction of the EGRB 
spectrum 26] with the recent Fermi data, and made 
quantitative predictions for the next four years of the 
Fermi mission. In Ref. [26| j . we constructed a model 
of the blazar gamma-ray luminosity function (GLF), 
taking into account the blazar SED sequence and the 
LDDE luminosity function inferred from X-ray obser- 
vations of AGNs. The GLF model parameters are 
constrained by fitting to the observed flux and red- 
shift distribution of the EGRET blazars. With this 
model, we can predict the EGRB spectrum in a non- 
trivial way, rather than assuming simple power-law 
spectra. 

The contribution from non-blazar AGNs to the 
EGRB is also considered, using the nonthermal coro- 
nal electron model [28| . This model gives a natu- 
ral explanation for the observed cosmic MeV back- 
ground, and we examined whether the X-ray to GeV 
gamma-ray background radiation can be accounted 
for by the two population model including blazars and 
non-blazar AGNs. 

Our predicted blazar EGRB spectrum matched 
very well with the latest Fermi EGRB data. This 
indicates that blazars are the dominant population of 
the EGRB and that other components such as dark 
matter annihilation do not contribute to EGRB signif- 
icantly. This also indicates that the use of blazar SED 
sequence was valid. The predicted spectrum including 
blazars and non-blazar AGNs is also in good agree- 
ment with the observed data from X-ray to 100 GeV. 
Therefore, blazars and non-blazar AGNs are likely the 



primary sources of the cosmic X-ray and gamma-ray 
background radiation. Our model predicts that it is 
possible to account for 80 % of the observed EGRB 
flux above 100 MeV with the sum of the blazar and 
non-blazar AGN components, where blazars alone ac- 
count for about 45 %. The two components are com- 
parable at ~100 MeV, with blazars and non-blazar 
AGNs dominating at higher and lower energies, re- 
spectively. 

From the viewpoint of the EGRB spectrum from 
MeV to GeV, we argue that the dominant component 
contributing to the MeV background is non-blazar 
AGNs that are also responsible also for the CXB, 
rather than MeV blazars as proposed by Ref. [29j . 

We predicted the flux distribution of blazars by our 
model, and we found that 750 and 1200 blazars in 
all sky should be detected by Fermi, assuming a sen- 
sitivity limit of Fi im = 3 x 10~ 9 and 1 x 10~ 9 pho- 
tons cm~ 2 s _1 above 100 MeV corresponding to one 
and five year survey sensitivity. This number is sig- 
nificantly lower than the predictions by many of the 
previous studies. About 4-50 of non-blazar AGNs are 
expected to be detected by Fermi with 5-year sur- 
vey. Fermi should resolve almost all of the EGRB flux 
from blazars into discrete blazars, with percentages of 
>98% with 5-year survey. On the other hand, less 
than 0.1% of the EGRB flux from non-blazar AGNs 
can be resolved into discrete sources. Therefore, we 
have a clear prediction: Fermi will resolve almost all 
of the EGRB flux into discrete sources at photon en- 
ergies > 1 GeV where blazars are dominant, while a 
significant fraction of the EGRB flux will remain un- 
resolved in the low energy band of < 100 MeV where 
non-blazar AGNs have a significant contribution. This 
prediction can easily be tested, only with the source 
counts and the EGRB estimates by Fermi data. 

We also made predictions for future detections of 
high-redshift blazars. Since the GLF is uncertain at 
z > 3 due to lack of direct observational data, we con- 
strained the redshift evolution of our original AGN 
XLF using SDSS quasar data. With this assumption, 
we find that Fermi may find some blazars up to z ~ 6 
during the five year survey. Such high-redshift GeV 
blazars will offer a new probe of the formation histo- 
ries of early stars and galaxies via GeV spectral at- 
tenuation signatures caused by the high-redshift UV 
background radiation. 
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